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Abstract 

The Eady problem of baroclinic instability as applicable to quasi- 
geostrophic oceanic flows with zero internal PV gradients is revisited 
by introducing a mild slope and Ekman pumping on the lower bound- 
ary. The solution behaviour is determined by the isopycnal slope 
relative to either the bottom slope or the ratio of Ekman depth to 
horizontal wavenumber. Attention is paid to the physical interpreta- 
tion of the growing, decaying and stable disturbances, with emphasis 
on the intimate connection between the quasigeostrophic edge waves 
and Eady waves, and the role of the isopycnal slope for the stability 
properties as opposed to the bottom density gradient. The distur- 
bance structure is found to be strongly influenced by the boundary 
conditions. 

For a sloping bottom boundary, the growth rate is enhanced for 
the most unstable waves if the isopycnals tilt in the same direction as 
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the bottom, but in general non-standard boundary conditions tend to 
retard the growth of disturbances. In particular, the existence of the 
long- and short-wave cutoffs is found to be very sensitive to boundary 
conditions, both for the sloping topography and the Ekman pumping. 
It is suggested that any cutoffs for the growth rate in an Eady-like 
problem actually result from the chosen boundary conditions. How- 
ever, for a certain range of parameters, the maximum growth rate is 
comparable to that found in the original Eady problem, which may ex- 
plain the fair success enjoyed by recent eddy parameterizations basing 
their timescale on the Eady growth rate. 

Keywords: Eady problem, baroclinic instability, topography, Ek- 
man pumping, eddy parameterizations 



1 Introduction 



The importance of oceanic transports induced by mesosc ale eddies has been 
demo nstrated in the previous decade for passive tracers nDanabasoalu et al 



1994) as well as, tentatively, for ecological, reactive tracers (\OschMes and Garcon . 
1998T ). The formation of these mid-ocean eddy transports is fr equen t ly de- 
scribed with the baroclinic, quasigeostrophic stability problem o uEadv (1949) 
as a prototype. 

Less attention has been devoted to eddy transports in the coastal areas. 
As opposed to the original problem formulated by Eady with only horizon- 
tal boundaries, coastal areas typically have a sloping bottom and, possibly 
as a consequence, buoyant coastal currents may persist over long distances 
under apparently unstable conditions. The stability of such currents is a 
difficult problem that hith erto has defied unapproximated analysis; but see 
Poulin and Swaters (1999) for recent progress. Extensions of the tractable 
analysis due to Eady are pursued in this study as a source of guidance for 
understanding the instability and in particular, the consequent eddy fluxes 
associated with coastal currents. 

In an atmospheric set ting, the effect of a smal l botto m slope in the Eady 
problem was stu died bv \Blumsack and Gierasch ( 1972j ) as a prototype for 
Martian climate. Mechosk (|l980al lbl) extended the theory to deal with two 
sloping boundaries and app lied it to the terrestrial at m osphere, w i th an 



de SzoekA (|l975h 



oceanographic discussion in \Mechoso and Sintonl ([19811 b 
developed the bottom profile into a perturbation series and discussed the 
resulting hybrid instabilities. 
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These studies found substantially modified stability properties for changes 
of the bottom slope, and as demonstrated at the end of the present investiga- 
tion, it is possible to relax the necessary condition for instability altogether 
to make the flow unconditionally stable. For the eddy parameterizations 
presently used in numerical models of the ocean, the growth rate is deter- 
mined from the original Eady theory. Hence, such parameterisations do 
not appear to be immediately applicable to coastal situations with sloping 
boundaries. 

The modifications ca used by Ekm a n laye rs at the bottom boundary layer 



are not obvious either. \Holopainera (|f 9611 1 discusses the effect of Ekman 



pumping in a two-layer model on a /3-plane, noting that friction may broaden 
the unstable wave-band while s till retaining the short- and long- wave cutoffs. 



Williams and Robinsonx (| 19741 ) generalized the Eady problem by combining 



the effects of a varying stratification with Ekman pumping at one or two 
boundaries. These authors found a destabilizing effect on the short waves, 
i.e. there is no parameter range with absolute stability. 

The instability in the Eady problem is often discussed as resulting from 
the meridional density- or temperature-gradient at the boundary. While 
this argument is mathematically correct, the actual reason for the boundary 
gradient is the intersection of sloping interior isopycnals with a horizontal 
boundary. The intimate connection between the topographically modified 



Eady waves and the quasigeostrophic edge waves touched upon by whines 



(1970) is in the present context found to be important for the stability prop 
erties. 

In the present work the growth rate and structure of the disturbances in 
baroclinic instability of the Eady problem in the pre sence of topography and 
frictio n are studied in more detai l than previously ( cf . \Blumsack and Gierasch 
( 1972[ ). Williams and Robinson ( 1974[ )). emphasizing the importance of the 
isopycnal slope for the physical interpretation. The lower boundary is mod- 
ified to include a slope (Section EJ) and Ekman pumping (Section OJ) in order 
to clarify their effects on the oceanic stability problem. 

The structure and stability properties of the solutions are found to be 
substantially modified by the boundary conditions and are discussed in Sec- 
tion 0] with emphasis on the physical background. From the analysis the 
pure Eady problem emerges as a special case of a more general baroclinic 
instability problem. 
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Figure 1: The Eady problem as applied to a buoyant flow with an infinites- 
imal bottom slope. The isopycnals of the basic-state density field p(y, z) 
slope with a constant angle throughout the fluid interior. Rigid lids, one 
horizontal on the top and one sloping on the bottom at h(y) <C H, prevent 
vertical movements across them. The flow U(z) has constant shear, with 
U(z = 0) = 0. 



2 The Eady problem with boundary layer slope 



Th e configurati on characterizing the continuous stability problem formulated 
by lEady (J1949) with a sloping lower boundary added is depicted in Figure 
^ For the atmosphere, an analytical ap proach to the problem has been 
pioneered bv Wlumsack and Gierascm ( 19721 ). whose analysis is extended here. 
Dimensional forms of the equations are retained in order to keep close t rack o f 
the physics. The relevant scaling arguments can be found in e.g. \Gill (1982). 
whose notation we follow; a brief derivation of the governing equations is 
presented to assist the discussion in Section HI 

The conservation equation for pseudopotential vorticity q in an incom- 
pressible flow reads 



D 
Dt 



Q2 Q2 Q 

— - H H 

dx 2 dy 2 dz 



N 2 dz 



D 
Dt 



q = 0. 



(1) 



Here ^ = ^ + + v ^ is the geostrophic advection operator with U 



[— fpo)~ 1 dp/ dy as the basic-state flow in the positive x-direction, v = (/po) _1 ^p/ dx 
is a geostrophic velocity in the y-direction, / is the Coriolis frequency, g the 
gravitational acceleration, p = p + p(y, z) the density, p the pressure field 
and iV 2 = — is the buoyancy frequency. 
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For a small perturbation q' (q = q{y) + q'(x, z, t)) on a time-independent 

background PV distribution q the linearized stability problem is formulated 

as _ 

D / dq 

with the perturbation velocity v. 

To keep the formalism straightforward, only an interior solution at a 
distance from the boundaries in the y direction is considered. The interior 
region is wide enough to permit the growth of waves and to justify the absence 
of Kelvin-wave-like boundary dynamics, but sufficiently narrow to justify the 
approximation of a small relative change in bottom depth: h = h(y) ^ H . 

The mathematical problem posed by Equation (J2J) becomes tractable if 
there are no gradients in the basic-state PV field. In this case dq/dy = 0, 
U = Uqz and the conservation law of Equation (J2J) becomes 

^ = 0, (3) 

i.e., q' is constant following a parcel. 

On the upper boundary, the condition for the vertical velocity is w\ z= h = 
0, whereas at the lower boundary impermeability of the boundary yields 
w\ z= h = vdh/dy. For small deviations from the background state and adia- 
batic flow, the boundary condition is obtained from the buoyancy equation 
Dp/Dt = 0. With the geostrophic advection operator and the hydrostatic 
approximation || = —pg, we arrive at the equation to be satisfied on both 
horizontal boundaries: 

^W£| = 0. (4) 

Equation may be recast in a more illuminating form by defining the 
isopycnal slope as 

a = -{dp/dy){dp/dz)- 1 = fN^ = iRr 1 *, (5) 

where dU/dz has been evaluated geostrophically as g(pf)^ 1 ^ and Ri is the 
geostrophic Richardson number. 

From Equation (J3J) and an application of Equation (@J on both boundaries, 
a linear stability problem for a perturbation (p = p Q (y,z) +p'(x,z,t)) is 
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obtained: 



dt + ~T Z dx~) \Jh? + lh [jPlh)) ~ ^ ' 

f d 2 p' d 2 p' dp' 

+ — — a Tr- = z = H (7) 



A^ 2 dzdt dzdx dx 
dh dp' ^ f d 2 p' ^ d 2 p' dp' ^ ^ . . 

dy dx N 2 dzdt dxdz dx 
The lower boundary condition should be evaluat ed at z = h, but for small 



Rossby numbers it can be linearized around z = ^Blumsack and Gierasch . 



19721 ); otherwise the baroclinic and barotropic parts are not separable. In 



dimensional form, the condition is that h be small compared to H. At this 
level of approximation, the bottom slope does not play a role for the hori- 
zontal momentum balance, and thus Kelvin-like boundary waves are filtered 
out from the problem. 

Taking A^ 2 constant, the variables are separated by writing 

p' = <f)(z)e ik{x - ct \ (9) 

which results in eigenvalue equations for the equation set (pj|)-([5]): 

<" -'>(£§?-*•*) = (10 » 

/ d6(0) (dh \ N 
We look for solutions to Equation (fT0|) of the form 

kNz k.Nz 

<P = C ie -— + C 2 e~ (13) 



Non-trivial solutions for the system of equations ([lip- ([T2]) must be found 
when the determinant of the associated coefficient matrix is zero. This con- 
dition yields a somewhat complicated characteristic polynomial, which can 
be simplified by defining a number of non-dimensional variables. 

The Rossby height Hr = f(kN)' 1 can be used to define the height ratio 
r/j = HH^ 1 = HkNf^ 1 , which is also the non-dimensional wave number. 
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The inverse of the buoyancy frequency times the isopycnal slop e, (vV| cr | ) 1 , 
appears as a natural time scale for the problem (ci. \Rhines ( 1970f )) and should 



be much larger than / _1 for the approximations leading to Equation (Q) to be 
valid. With this time scale, the frequency uj = ck can be non-dimensionalised 

as 

U ~ No ~ fRiV 2 ~ H R dU/dz' [ j 

One possible physical interpretation of uj is as the ratio between the phase 
speed and the background current velocity at one Rossby height above the 
bottom. 

A further notational simplification is obtained by defining the ratio be- 
tween the bottom and isopycnal slopes as A = |^cr _1 , which leads to the 
following characteristic equation: 

(1-A) ((l-e 2 ^) + (l + e 2rh H) + 
+ (A(l + e 2rh ) + (1 - e 2rh )r h ) uj - (1 - e 2rh )uj 2 = (15) 



with the solution for uj is obtained as 

1 

2(e 2 ^ - 1) 



& = «, J n ( - A ( e2rh + 1) + ^{e lrh - 1) ± 



± a/ (2 - A + r h y + e 4 ^ (-2 + A + r h y + 2e 2 ^ (-4 + A (4 + A) - r h 2 ) .(16) 



2.1 Solution properties 

The standard Eady problem is recovered by setting A = 0. It is readily found 
that in this case imaginary values are possible only for < < 2.39936, i.e. 
the upper and lower boun daries must not be more than about 2.4 Rossby 



heights apart (\Eadii Il949h . 

With a sloping bottom boundary, the imaginary part of uj for a range of 
interesting A- values is given in Figure 12(a), where unstable solutions only 
are possible within the contoured area. For opposing isopycnal and bottom 
slopes, there is a strong stabilisation of long waves, while the short-wave cut- 
off moves to even shorter waves. A narrow unstable region remains, however. 

For slopes in the same direction, the maximal growth rate initially in- 
creases as A grows, before decreasing as A — > +1 (the maximum of uj is 
0.366904 for r h = 1.13777, A = 0.398331, as opposed to 0.309817 in the 
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12345 012345 
(a) H/H * (b) 

Figure 2: (a) Non-dimensional growth rate = ckN^fdU/dz) -1 as a 

function of the height ratio = HH^ 1 = HkNf -1 and the slope parameter 
A = |^er _1 . Contours 0, 0.1, 0.2 and 0.3 are drawn, (b) Non-dimensional 
frequency (9ft(cJ+)) as a function of and A. Contours increase from zero in 
the upper left hand corner with 0.2 intervals. 



Eady case when = 1.60609). For A > 1, i.e. when the isopycnal slope 
is smaller than bottom slope, no unstable solutions are to be found, a result 
that is formally demonstrated in Section 14.21 

The real part of to in Figure 12(b) also shows intriguing characteristics 
and seems to have escaped discussion . (Fig ure ^a) is essentially similar 
to Figure 2 of Blumsack and Gieraschl (|l972l ) but these authors appear to 
have neglected the cont ribution from th e discriminant to the real part of 
the solution. Note that de Szoeke\ ( 1975f ) gives a slightly different stability 
diagram, apparently because of his bounded y-domain.) It is readily seen 
that in the unstable wedge, the waves are only weakly dispersive since to 
increases in an approximately linear fashion towards oc k, i.e. c ~ const. 
However, outside the unstable wedge and disregarding the very longest waves, 
the waves are dispersive, implying that wave packages with different wave 
numbers would propagate away from each other. 

The different propagation characteristics are more understandable in the 
light of the structure of the disturbances (Figure |3J). It is recognised that 
the dispersion properties are determined by the boundary at which the dis- 
turbance is trapped; on the lower boundary, Cj oc A, whereas on the upper 
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boundary Co oc r h . 

The wave properties and the resulting instabil ity may be i nterpre ted as 
the interaction between two boundary waves (e.g. \BrethertoA (|l966bh . For 
uj + and opposing isopycnal and topographic slopes (A < 0), the long waves 
are confined to the lower boundary. For increasing wavenumbers, the upper 
boundary comes into play, and the perturbations on both boundaries start to 
interact and grow in time, because the shear maintains their phase difference 
while the mutually induced velocity field tends to diminish it. For even larger 
wavenumbers, the contribution from the lower boundary vanishes and what 
remains is a stable wave on the upper boundary. 

The most unstable solution differs somewhat in character from the Eady 
wave, since it is more confined to the upper boundary. This is caused by the 
sloping bottom boundary amplifying the vorticity contributions for small per- 
turbations. The result is that in the interaction regime, larger perturbations 
on the upper boundary are permitted, resulting in an overall amplification 
of the instability. 

For cD_, the situation is reversed: long waves are found on the upper 
boundary, short ones on the lower boundary, and in the interaction zone the 
upper boundary perturbation is "upstream" of the lower boundary, which 
gives rise to a destructive interaction. 

Moreover, there is a narrow regime of negative frequencies for 0.4 < A < 
0.6 and 0.4 < < 0.6. The possible relevance of this regime to a more 
realistic problem has not been discussed, but it appears as though the edge- 
wave phase speed exceeds the Eady-wave phase speed in this regime, resulting 
in wave progression with the coast on the right. 



3 Eady problem with Ekman pumping on the 
lower boundary 

The effects of a viscous boundary layer on the interior balance may be repre- 
sented by diagnosing the vertical velocity at the top of the bottom boundary 
layer using the steady-state Ekman transport, 

r v\ 1/2 fdv _ duA j_ / v\ 1/2 / ay gy \ 

WE ~\2f) \d~x dy)~pf\2f) ydx^ + W) 1 ( } 
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Figure 3: The structure within a constant, i.e., p'Cf 1 , of growing (a) and 
decaying (b) perturbations in the (r h , A) -space. The vertical axis is scaled 
with the domain height, z = zH" 1 , and the horizontal is scaled with the 
Rossby radius, x = xL~j^ = xf(NH)~ 1 with a total extent of Arc. Outside 
the region of instability, the waves are confined to either boundary. Mean 
shear is to the right (positive x direction) in the figure 



10 



which in the absence of y- variations can be written as we = jj f^jj gjr- 
Hence, the lower boundary condition of Equation (JBJ) can be reformulated as 

v \ 1 d 2 p' f d 2 p' d 2 p' dp' „ 



2/ / dx 2 N 2 dzdt dxdz dx 

The solution procedure of Section |2 can now be repeated, resulting in a 
partly imaginary characteristic polynomial 



e 



2rh (2-iA„-2a>) (-l + r ft -cD) + (l + r ft -cD) (2 - i A„ + 2a)) = 0,(19) 



where A v = ko~ l ^Jvf~ l = E v o~ x . E v = k\/ vf~ x is an Ekman number 
based on the horizontal wave number rather than the depth scale; note that 
E v is also a slope-like parameter in that it is the ratio of a vertical to a 
horizontal lengt h. A„ is related to the ratio be tween the Ekman and Rossby- 



Kibel numbers ^Williams and Robinsom Il974l cf. their definition of Q) 



The solution of Equation (fTU|) is 

Cj = (4 (e 2r » - I))' 1 (iA v (e 2rh + l) - 2r h [e 2 ^ - l) ± 
i^-A, 2 ^" + l) 2 + 4iA„ {e 2rh - 1) (2 + e 2r * (r h - 2) + rj + 

+4 (e 2 ^ - 1) (e 2 ^(^ - 2) 2 - (r h + 2) 2 )^) ' J (20) 

The associated stability properties are shown in Figure 0] Note that 
since a change in the sign of the isopycnal slope affects both A and the 
dimensional version of u (Eq. ((HI)), the A u < part of the uj + solution is 
actually decaying, and vice versa for cu_; the positive and negative solutions 
switch roles when A u changes sign and hence only the A v > region is 
shown. 

The stability characteristics of u + in Figure 0^a) can again be understood 
from the structure of the perturbations, as presented in Figure E(a). For 
the short waves that are stable in the Eady problem (77, > 2.4), there is a 
change in the character of the solution with changes in A u . The waning Eady 
solution, which is bottom-trapped in the short-wave range, is barotropised 
and becomes strongly decaying by increased Ekman pumping. The damping 
increases exponentially in the long-wave limit — > 0+. 
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(a) H ^ (b) W 

Figure 4: Growth rate for (a) the growing solution 3[d;+] and (b) the de- 
caying solution to the Eady-Ekman problem, as a function of the 
non-dimensional friction A v and boundary separation = HH^ . The de- 
caying solution is obtained with negative A„ due to the symmetry of square 
roots of complex numbers. 



The growing solution is top-trapped and destabilized for short wave- 
lengths due to baroclinisation. Moreover, in the range > 2.4, the maxi- 
mum growth rate is always found for non-zero values of A u , since the Eady 
solution A v = is stable. 

The stability properties also manifest themselves in the phase speed of the 
waves as an increase (decrease) for decaying (growing) long waves, cf. Fig- 
ure ini 

4 Discussion 

In the previous sections the effect of the lower boundary condition on the 
phase speed of baroclinic disturbances and their vertical trapping was found 
to be substantial. Keeping in mind the basic mechanism of instability in the 
Eady problem, we propose that any non-standard boundary condition that 
modifies the boundary-wave phase speed will affect the stability properties in 
a rather fundamental way. Some possible consequences of varying boundary 
conditions and their interpretations are discussed in the following subsections. 
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Figure 5: Structure of (a) positive (u) + ) ( (r h = 1.6, A„ = 5), (4.0,2.0)), 
and (b) negative solutions ((1.6,1) and (4.0,2.0)) in the Eady-Ekman 
problem. Note the change in x-axis length for rightmost figures. 




(a) m (b) m 

Figure 6: Real part of Co for (a) the growing and (b) the decaying disturbance 
in the Eady-Ekman problem. 
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4.1 The edge wave — Eady wave duality 

The phase speeds of the solutions to the bottom-slope problem when the dis- 
turbances are trapped to either boundary may be approximated from Equa- 
tion (JTfij) as 

lim uj + = — 1 (21) 
lim w_ = 1 - A, (22) 

which yields the following dimensional results: 



Cupper 



(dU N 2 dh\ rI , dh N 

<*—=(&- = a,)*- < 24 > 

This result for non-interacting boundary waves could also have been obtained 
by inserting the exponential solution of Equation (JT3j) separately into Equa- 
tion (|T2|) or Equation (|TTjl and solving for c, yielding a wave solution for one 
boundary only. 

For the topographically modified Eady problem with positive shear, the 
lower boundary phase speed is reduced by a boundary with a positive slope, 
because the slope generates a quasigeostrophic edge wave that tends to travel 
in the Kelvin-wave direction, opposite to the Eady wave and therefore decel- 
erating it in the case of positive shear, but accelerating it if A < 0. Indeed, 
even negative phase speeds are possible on the lower boundary. 

Within this approximation, there is a duality between Eady waves with 
sloping isopycnals on a flat bottom and e dge waves with horizontal isopycnals 
on a sloping bottom (cf. \Rhines ( 1970h ). The same phase speed as with 



dh/dy = 0, a = (T f° r Qower is obtained with o = if dh/dy = —cr . The 
duality is illustrated in Figure For o = 0, there are no unstable solutions. 

Note also, that the phase speed on the lower boundary in the absence 
of an isopycnal slope, ci ower = — f^y, is frequently referred in the literature 
as the Nof speed, being the speed at which eddies progress on a sloping 



bound ary. An explicit connection to the edge wave phase speed of \Rhines 
is appropriate in this context. 
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(a) (h) V> 

Figure 7: Duality between stable Eady waves with a flat bottom but sloping 
isopycnals (a) and edge waves with flat isopycnals but sloping bottom (b): the 
waves on the boundary have the same phase speed, only the angle between 
the boundary and the isopycnals matters 



4.2 Surface temperature gradient vs. isopycnal slope; 
the necessary condition for instability 

It is customary to decompose the quasigeostrophic dynamics into contri- 
butions from the inter i or PV a nd from the boundary tempera ture gradi- 
ent (e.g. \Hoskins eta}] (jl985h). \ffnlnpaine.n and Kaumla\ fll99lh have crit- 
icized this decomposition by arguing that the quantities depend on each 
other, and when evaluated from observational data there is a high degre e 
of cancellation in the inv e rted fi elds over large scales (see also \Davis (1993); 
Holopainen and Kaurola ( 19931 )). We briefly revisit this issue from a stability 
perspective. 

The vertical divergence d/dz of the boundary condition of Equation P|) 
may be reformulated as 



Tft + U[z) lk 



d f 2 dp 
dz N 2 dz 



-—( f w ) - y—— ^ dp 

dz dy dz N 2 dz ' 



(25) 



which, with w = vdh/dy, yields 
D d ( pdp 



Btdz \N 2 dz 



+ fph = 



0,H. 



(26) 



The range of validity of the Poisson equation (JHJ) may now be extended 
to the full vertical domain, if the boundary conditions are generalised to a 
flux, which is constant in the interior and diverges on the boundaries, so that 

the y-term of the material derivative becomes v-^-^H(z) (^jj^ + f 2 phj = 
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vJt-5(z) (jJrQ- + Pph). Here H(z) is the Heaviside step function and its 



\N 2 dz 

vertical derivative, the 5-function, ensures that the boundary condition only 
appears exactly at the boundary. Thus the boundary effects may be incor- 
porated in Equation ((TJ) as 



D 

Wt q 



D 

m 



D 



S(z] 



f fdp 
\N 2 dz 



+ f 2 ph 



d 2 p 
— - + 

dx 2 

-5{z 



d_ 

dz 



f_dp 
N 2 dz 



} N 2 dz J 



(27) 



Brethertonx (|1966al) was the first author to suggest a form similar to Equation 
(J2ZJ), but without considering it as the divergence of a flux. 

The form of PV in Equation ()27|) is c onvenient when app l ying i ntegral 
theorems, such as the stability theorem oi lCharnev and Stern ( 1962f ). This 
theorem requires a change in sign of dqj dy for instability. Since there are no 
interior PV gradients in the Eady problem, the ^-gradients in t he boundary 
PV co ntributions can be recognized as the source of instability \Bretherton . 
1966al ): compare Equation (|2*7j) with Equation (0). 

It is now possible to pose the question whether the source of instability 
at the lower boundary can be eliminated. This proves to be the case when 
both contributions at z = add to zero in Equation (j2*7|) : -^s§^ = ~f 2 P^> 



i.e. when 



dz dp 
dp dy 



dh 
dy' 



(28) 



This shows that the removal of the necessary condition for instability is 
possible when the the bottom slope is equal or larger than the isovycnal slope. 
This p ossibility for stabilisation is also remarked on bv \Charnev and Flierl 
(1981). The appropriate physical explanation, however, is found by recog- 



nizing that the boundary density gradient, which is often referred to as the 
source of instability, is a side-effect of the isopycnal surfaces intersecting the 
horizontal boundary. 

In the stabilized case, the isopycnals intersect only the top boundary; they 
do not feel the vorticity effect of the bottom because they are "parallel" to it. 
Hence, by removing the lower boundary PV source, the necessary condition 
for instability is also removed, an d the flow is unco nditionally stable, as seen 
in the growth rates in Figure 121 $Mechoso\ (|1980a[ ) explained this removal of 
instability using heat-transport arguments, but since the quasi-geostrophic 
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dynamics are determined by vertical velocity divergence, the heat fluxes are 
a consequence rather than a cause of the stabilisation.) 

This result, strictly valid under the narrow limitations of quasi-geostrophy, 
makes a lot of physical sense and is therefore expected to be relevant also for 
non-quasi-geostrophic problems. Indeed, a parameter similar to A a ppears 



independently in several contexts in the literature, e.g. \Swaterd (|1991[ ) found 
a similar condition for a non-quasi-geostrophic frontal model. However, other 
types of instabilities that may emerge under different approximations cannot 
be ruled out. 

For opposing isopycnal and bottom slopes, the necessary condition for 
instability is fulfilled, but nevertheless considerable stabilisation of in partic- 
ular long waves takes place already for very small bottom slopes, cf. Figure 
® 

It is useful to note that the source of vorticity on the boundaries is due to 
their impermeability, implying w — 0. The divergence of vertical velocity on 
the boundaries is the source of the 5-function vorticity sheet, which makes the 
angle between the isopycnals and the bottom the relevant parameter for the 
phase propagation. The validity of this argument is born out by examining 
Equation (0J), which may be reformulated as 

p/ 2 ^ = --f#4^V (29) 



dz m\dzN 2 dz / 

For the advecting ?/-part of D/Dt, Equation (@J) yields w = va, i.e. the flow 
of v against the isopycnals causes a vertical velocity that diverges only when 
a changes or a boundary is encountered, as illustrated in Figure |H1 This 
divergence introduces a change of the interior PV, as indicated by (J2HJ). For 
the linearized problem, the isopyc nal surfaces of t he ba sic state are imper- 
meable. This isentropic upgliding ( Hoskins et al . 1985| ) is exactly canceled 



at the boundaries in the Eady problem by the boundary- induced velocity 
("boundary temperature advection") to make the net w zero. 

From a physical standpoint, it thus makes a lot of sense to divide the 
quasi-geostrophic flow into contributions from the horizontal vorticity (Vjjp), 
and the modifications caused by divergence of vertical velocity (^^2 + 5- 
functions), giving rise to the familiar phenomenon of vortex stretching. The 
natural corollary of this argum ent is the second flow partition suggested by 



Holopainen and Kaurolca ((1991). 



We further remark that it is customary to draw para llels betwe e n the 
stabilisation caused by a bottom slope and the (3 effect ( Junaclaus . 1999; 
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Figure 8: Illustration of the divergence of vertical velocity in the Eady prob- 
lem with a sloping bottom. A horizontal flow against isopycnals is converted 
into vertical velocity, which diverges where isopycnals intersect the bound- 
aries, generating vorticity by vortex stretching. The bottom slope creates a 
divergence of its own. 



Spall and Chapman . 1998| ). There is, however, a clear difference between the 



planetary vorticity gradient and the effect of topography. The bottom slope 
appears in the purely baroclinic problem as a 5-function vorticity contribu- 
tion at the boundary, whereas the /3-effect is homogeneous over the whole 
depth. Therefore, the effect of the slope on the baroclinic instability is fun- 
damentally different from the /3-effect, whereas for the barotropic instability 
the effect is the same, but only within the linearisation permitted by h -C H. 
Furthermore, as Charnev and Flierll ( 1981 ) point out, the variation of N with 
depth in the ocean also makes the vorticity dynamics (and hence, the sta- 
bility properties) deviate from those of the simple model based on constant 
stratification. 



4.3 The sensitivity of eddy parameterisations 

The recent decade has seen a rise of interest in parameterising the effect of 
eddies in coarse (viz. non-eddy-resolving) ocean models. Already a rather 
coarse AGCM is capable of resolving the atmospheric eddies, whereas an 
eddy-resolving OGCM is a major computational challenge that will be unaf- 
fordable for routine research for many year s to come. 

Building on the atmospheric studies of Green ( 1970j ) and Stone 



18 



Visbeck et al\ ()1997f ) based their eddy pa ramet e risatio n on the baroclinic in- 
stability growth rate in the problem of Eadv (1949). Subsequently, many 
studies have found that a parameterization of eddy fluxes based on the Eady 
growth rate can be brought to reasonable harmony with ed dy-resolving 3 



dimensional model results far from the lateral boundaries ( Visbeck et al. 



1997c \Jones and Marshall. Il997t \Spall and Chapmanl 11998), though only in 



Gille and Davis 



an en semble-mean sense and with significant error bars 
1998h . 

In light of the results from Sections EJandElwe may now inquire about the 
sensitivity of the parameterisation schemes referred to above on the details of 
the lower boundary, with emphasis on variations ty pically foun d in coastal- 
flow situations. Based on dimensional arguments (cf. \Held\ (|l99Sj) for a recent 



overview), an eddy- induced diffusivity is found to have a timescale T and a 
length scale L associated with it such that 



K OC — — OC 

T 



f 



-r 



(30) 



The in verse of the instability growth rate (fyBA , as used by 1 Visbeck et al. 



( 1997t ) ) is a natural timescale in the eddy diffusion process, but the length 
scale is mor e ambiguous. L may be determined by internal dynamics (cf. 
Stone\ (|l972| )). or by th e external geo metry of the problem (e.g. the width of 
the baroclinic zone, cf. Green ( 1970| )). Hence, 



K int ^ ji 



w max 



/i,max 



or n ext oc ^max^ext- 



(31) 



Since a variation in A or A„ is expected to affect the internal length scales, 
the inverse of the fastest-growing wavenumber has been taken to represent 
the internal length scale for Note that L ext is independent of A and 

A u , which implies a difference between and ft ex t- Furthermore, because 
of the dependence of cDmax on A, the diffusivity and hence the eddy flux will 
be sensitive to A. 

The dependence of k on the bottom slope and Ekman pumping (A) is de- 
picted in FigureEta). The largest deviations from the standard Eady problem 
are encountered when the eddy length scale is set by the fastest growing wave; 
in the diagram k is vary i ng wit h a factor of 10 in the range —0.5 < A < 0.5. 
Blumsack and Gierasch ( 1972h calculated the correlation v'T' as a function 
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K K 




Figure 9: (a) the eddy diffusivity K- m ^ (solid) and n ex ^ (dashed) for the slop- 
ing bottom (b) influence of the bottom Ekman layer on the eddy diffusivity 
based on the internal and external length scales, as in (a). The diffusivities 
corresponding to the original Eady problem are normalised to unity. 



of A and obtained a result which shows a pronounced similarity to our Kj n ^, 
hereby providing support for our simple scaling argument. 

However, if the diffusivity is governed by a fixed external length scale, k is 
proportional to uJm&x and is therefore less sensitive to A. In case of bottom 
Ekman pumping, the sensitivity of k to A is smaller than in the presence of 
a sloping bottom, with k varying only by a factor of two for < A„ < 1. 

Hence, when the deviations from the Eady problem are not large, it ap- 
pears that the bottom slope exerts a strong influence on the eddy diffusivity 
when the length scale is determined by internal dynamics. Furthermore, the 
possibility for complete stabilisation in the case of a sloping bottom (imply- 
ing zero diffusivity), and the destabilisation of short waves in the Ekman 
layer case call for a development of new eddy parameterisations for coastal 
waters. 

5 Conclusions 

The baroclinic instability mechanism requires the interaction between two 
propagating PV disturbances. In the Eady problem these disturbances re- 
sult from the divergence of vertical velocity at the solid boundaries, and are 
therefore boundary-trapped. The relation between the isopycnal slope and 
the bottom slope (or ratio of Ekman depth to horizontal wavelength) is the 
single most important parameter in this problem. 

The propagation speed of the disturbances influences the stability prop- 
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erties of the problem. An impermeable boundary which is not horizontal 
leads to changes in the phase speeds of the disturbances. In particular, the 
existence of short- and long-wave cutoffs appears to be especially sensitive 
to the prescribed boundary conditions via the modifications these give rise 
to in the boundary phase speeds. 

Therefore it seems that an accurate representation of the Eady instability 
growth rate in an eddy parameterization would require quite detailed knowl- 
edge of the phase speeds of free Eady waves on horizontal boundaries. These 
phase speeds, on the other hand, are sensitive to the details of the boundary 
even in the quasi-geostrophic formulation. One may expect this sensitivity 
to be carried over to mildly non-quasi-geostrophic dynamics. 

The importance of the trapping boundary for disturbances is emphasized 
in this study. In the Ekman-pumped Eady problem, all amplifying modes 
are top-trapped and inclined against the shear in the Eady amplification 
sense. The effect of Ekman pumping on the long waves is to reduce the 
growth rate, but the same structural change causes a mild destabilisation 
for short waves. The pumping translates the vorticity associated with the 
top disturbance into a divergence in the bottom layer, thereby introducing 
a small disturbance at the bottom. The top disturbance and its "slave" on 
the bottom have a phase shift, which is responsible for the destabilisation. 
Remarkably enough, a single disturbance in combination with a suitable 
boundary condition suffices for the generation of instability. 

The fairly large changes in the disturbance structure for bottom slope 
or Ekman pumping are probably a useful guide when analysing observations 
and model results, and since these effects are almost inevitably present in 
laboratory investigations, numerical models and observations, the classical 
Eady problem may not be the best guide when interpreting experimental 
results. 

The solutions derived in this study indicate that idealized numerical eddy- 
transport experiments without bottom slope or Ekman pumping may not 
yield universally valid parameterizations. As has been shown here, there 
are regions of parameter space where the stability properties of the quasi- 
geostrophic problem differ markedly from those found in the basic Eady 
problem. Most importantly, there is a possibility for the removal of the 
necessary condition of instability when the isopycnals are aligned with the 
bottom slope. For larger slopes, Kelvin-wave dynamics become important 
and the stability properties may be expected to differ from those studied 
here. 
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However, for small slopes and small Ekman numbers, the maximum 
growth rates do not differ very much from those found in the unmodified 
Eady problem. This probably explains the reasonably high level of success 
enjoyed by eddy parameterizations based on the Eady growth rate. 

The vertical-velocity and isopycnal-slope perspective applied throughout 
th is study serves as an interesting complement to the "IPV thinking" due 
to Hoskins et al. The straightforwardness with which the stabilizing 

influence of topography can be understood within the present conceptual 
framework is a good argument for considering the interior as well as the 
boundary contributions to vortex stretching in terms of divergence of verti- 
cal velocity, as opposed to the less intuitively evident ^-sheets of PV. Since 
boundaries have such a strong influence on the evolution of baroclinic pro- 
cesses, a robust specification of the boundary condition would be important 
for inversion purposes. In addition, the specification of a relevant boundary 
density distribution in the real atmosphere or ocean may be more difficult 
than the specification of large-scale isopycnal slope for use in the inversion 
process. 

The parameterisation derived by scaling arguments shows a strong de- 
crease of eddy diffusivity with an increasing bottom slope when the mixing 
process is governed by an internal length scale. Such a decrease is consis- 
tent with the persistence of buoyant coastal currents. However, the quasi- 
geostrophic theory does not appear to fully explain the observed persistence, 
since even with opposing isopycnal and bottom slopes, there is a narrow 
range of unstable wavenumbers. This unstable range would in principle be 
sufficient to cause a disintegration of such currents. But in fact, the unstable 
range is so narrow that a small modification of the boundary wave propaga- 
tion properties may be capable of stabilizing the flow in toto, thus bringing 
the theory into harmony with observations. Such a modification is, however, 
probably beyond of the scope of quasi-geostrophic analysis. 
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